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Report  Title 

Multi-linear  Stress  Strain  and  Closed-form  Moment  Curvature  Response  of  Epoxy  Resins 

ABSTRACT 

A  simplified  multi-linear  stress  strain  approach  has  been  used  to  obtain  the  closed  form  nonlinear  moment  curvature 
response  for  epoxy  resin  materials.  The  model  consists  of  constant  plastic  flow  in  tension  and  constant  yield  in 
compression.  The  multi-linear  stress  strain  model  is  described  by  two  main  parameters  in  addition  to  four 
non-dimensional  tensile  and  five  non-dimensional  compressive  parameters.  The  main  parameters  are  modulus  of 
elasticity  in  tension  and  strain  at  the  proportional  elastic  limit  point  in  tension.  The  nine  non-dimensional  parameters 
are  strain  at  the  ultimate  tensile  stress,  maximum  strain,  post  elastic  proportionality  stiffness,  and  post  peak  strength  in 
the  tension  model  and  strain  at  the  proportionality  elastic  limit,  strain  at  yield  strength  point,  maximum  strain,  initial 
elastic  stiffness  and  post  elastic  proportionality  stiffness  in  the  compression  model.  Explicit  expressions  are  derived  for 
the  stress-strain  behavior  of  the  polymers.  Closed  form  equations  for  moment  curvature  relationship  are  presented.  The 
results  of  tension,  compression,  and  bending  tests  using  digital  image  correlation  technique  are  presented.  Load 
deflection  response  of  flexural  three  point  bending  (3PB)  samples  could  be  predicted  using  the  moment  curvature 
equations,  crack  localization  rules,  and  fundamental  static  equations.  The  proposed  nonlinear  moment  curvature  shows 
good  predictions  when  compared  to  experimental  results. 
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Abstract 

A  simplified  multi-linear  stress  strain  approach  has  been  used  to  obtain  the  closed  form 
nonlinear  moment  curvature  response  for  epoxy  resin  materials.  The  model  consists  of 
constant  plastic  flow  in  tension  and  constant  yield  in  compression.  The  multi-linear  stress 
strain  model  is  described  by  two  main  parameters  in  addition  to  four  non-dimensional 
tensile  and  five  non-dimensional  compressive  parameters.  The  main  parameters  are 
modulus  of  elasticity  in  tension  and  strain  at  the  proportional  elastic  limit  point  in 
tension.  The  nine  non-dimensional  parameters  are  strain  at  the  ultimate  tensile  stress, 
maximum  strain,  post  elastic  proportionality  stiffness,  and  post  peak  strength  in  the 
tension  model  and  strain  at  the  proportionality  elastic  limit,  strain  at  yield  strength  point, 
maximum  strain,  initial  elastic  stiffness  and  post  elastic  proportionality  stiffness  in  the 
compression  model.  Explicit  expressions  are  derived  for  the  stress-strain  behavior  of  the 
polymers.  Closed  form  equations  for  moment  curvature  relationship  are  presented.  The 
results  of  tension,  compression,  and  bending  tests  using  digital  image  correlation 
technique  are  presented.  Load  deflection  response  of  flexural  three  point  bending  (3PB) 
samples  could  be  predicted  using  the  moment  curvature  equations,  crack  localization 
rules,  and  fundamental  static  equations.  The  proposed  nonlinear  moment  curvature  shows 
good  predictions  when  compared  to  experimental  results. 

Keywords  epoxy  resins,  stress  strain  relation,  moment  curvature,  material  properties, 
plastic  flow,  nonlinear  behavior,  flexural  response,  load  deflection,  three  point  bending 

Introduction 

Epoxy  resins  are  one  of  the  frequent  matrix  materials  in  fiber  composites.  Mechanical 
properties  (stress  strain  relationship)  and  progressive  failure  is  still  a  challenge  for 
researchers.  Difficulty  of  a  constitutive  law  in  polymer  matrix  composites  is  mainly  due 
to  the  characterization  of  polymer  mechanical  behavior.  The  hydrostatic  component  of 
stress  has  a  significant  effect  on  the  load  deformation  response  of  resins  even  at  low 
levels  of  stress.  Hydrostatic  stresses  are  known  to  affect  the  yield  stress  and  nonlinear 
response  of  polymers;  the  absolute  value  of  the  yield  stress  in  compression  is  higher  than 
the  ultimate  tensile  stress.  In  order  to  develop  a  general  model  for  polymer  composite 
materials,  the  behavior  of  polymer  resins  under  different  types  of  loading  has  to  be 
understood.  Wineman  and  Rajagopal9  (2000)  used  a  viscoplasticity  model  to  capture  the 
polymer  behavior.  Zhang  and  Moore10  (1997)  used  the  Bodner-Partom  internal  state 
variable  model  originally  developed  for  metals  to  obtain  the  nonlinear  uniaxial  tensile 
response  of  polyethylene.  By  modifying  the  definitions  of  the  effective  stress  and 
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effective  inelastic  strain  rate  in  the  Drucker-Prager  yield  criteria,  Li  and  Pan5  (1990), 
Chang  and  Pan1  (1997),  and  Hsu  et  al.3(1999)  developed  a  viscoplasticity  approach  for 
the  constitutive  law  of  polymer  materials.  Gilat  et  al.  (2007)  used  an  internal  state 
variable  model  to  modify  the  Bodner  model  to  capture  the  effects  of  hydrostatic  stresses 
on  the  response.  In  their  approach,  a  single  unified  strain  variable  is  defined  to  represent 
all  inelastic  strains.  Jordan  et  al.4  (2008)  modified  the  original  Mulliken-Boyce7  model 
(2006)  for  one  dimension  to  capture  the  compressive  mechanical  properties  of  polymer 
composites.  The  original  model  is  a  three  dimensional  strain  rate  and  temperature 
dependent  model  for  thermoplastic  polymers.  The  majority  of  the  parameters  were 
determined  by  fitting  the  model  to  experimental  compressive  data.  In  this  study,  the 
flexural  behavior  of  a  beam  is  investigated  in  an  attempt  to  establish  a  relationship 
between  the  tensile  and  compressive  stress  strain  curves  in  one  side  and  moment 
curvature  response  of  epoxy  resin  material  in  the  other  side.  In  order  to  correlate  tension, 
compression  stress  strain  curves  and  flexural  data,  a  closed  form  solution  has  been 
developed  to  obtain  moment  curvature  response.  The  load  deflection  response  for 
nonlinear  materials  under  determinate  static  conditions  has  been  developed.  Using 
inverse  analysis,  the  effect  of  stress  gradient  on  the  multi-linear  stress  strain  curve 
obtained  from  the  in-plane  uniaxial  tests  has  been  studied. 

Tension  and  Compression  Multi-linear  Stress  Strain  Curve 

The  multi-linear  stress  strain  curve  for  tension  and  compression  is  bilinear  up  to  the  peak 
stress.  Figure  1  shows  the  simplified  tension  and  compression  stress-strain  relationship  of 
epoxy  resin  materials.  The  tension  and  compression  curves  are  defined  uniquely  by  the 
parameters  E,  Spel,  Hu,  Hut,  oc,  co,  y,  f:L  jUco,  Hd,  and  Hue-  The  tensile  stress  at  the 
proportionality  elastic  point  (PEL)  is  related  empirically  to  the  stress  at  the  ultimate 
tensile  strength  (UTS)  point.  The  ascending  part  of  the  tension  and  compression  stress 
strain  diagrams  consist  of  two  linear  parts:  0  to  PEL,  and  PEL  to  UTS  in  tension  or  PEL 
to  compressive  yield  stress  (CYS)  in  compression.  The  curve  after  peak  strength  is 
idealized  as  horizontal  with  Gf  and  gcys  as  the  post  peak  sustained  stress  in  tension  and 
constant  yield  strength  in  compression  respectively.  The  constant  post  peak  tensile  stress 
level  co  shows  the  ability  of  the  model  to  represent  a  continuous  (co=  1)  or  discontinuous 
stress  response.  The  post  peak  response  in  tension  terminates  at  the  ultimate  tension  strain 
level  (sut  =  Hut  £pel ),  and  for  compression  it  ends  at  ultimate  compression  strain  level  (eUc 
=  Hue  Spel).  In  the  elastic  range,  the  resin  beam  in  bending  could  be  treated  as  a  bi¬ 
modulus  material  with  different  moduli  in  tension  and  compression.  The  tension  and 
compression  stress  strain  relationship  are  defined  as  shown  in  table  1. 


Table  1.  Multi- linear  stress  strain  curve 
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<jc ,  <jt,  sc,  and  s,  are  compression  and  tension  stresses  and  strains  respectively.  The  nine 
normalized  parameters  used  in  the  definition  of  the  constitutive  law  are  defined  by 
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Using  classical  beam  theory,  linear  distribution  of  strain  across  the  depth  is  assumed.  The 
stress  and  strain  distribution  across  a  section  of  a  beam  with  depth  h  and  width  b  by 
imposing  normalized  top  compressive  strain  in  different  cases  are  shown  in  Figure  2. 
Normalized  heights  of  compression  and  tension  sub-zones  with  respect  to  beam  depth  h 
are  shown  in  Table  2.  Tables  3  and  4  present  the  normalized  stress  at  the  vertices  of  the 
tension  and  compression  sub-zones  with  respect  to  tensile  stress  at  the  proportionality 
limit  point.  The  internal  force  in  each  compression  and  tension  sub-zone  of  nine  stress 
distribution  cases  is  calculated  from  the  stress  diagram  and  the  normalized  form  with 
respect  to  the  tension  force  at  the  proportionality  limit  point  ( bhEsPEL )  is  presented  in 
Tables  5  and  6.  The  centroid  of  the  stress  in  each  sub-zone  represents  the  line  of  action 
and  moment  arm  respect  to  the  neutral  axis  and  is  presented  in  the  normalized  form  in 
Tables  7  and  8. 


(a)  (b) 


Figure  1 .  (a)  constant  flow  in  tension;  (b)  constant  yield  in  compression 


(a)  rectangular  cross  section 


(b)  case  one 


(c)  case  two 


(e)  case  four 


(g)  case  six 


(d)  case  three 


(f)  case  five 


(h)  case  seven 


(i)  case  eight 


(])  case  nine 


Figure  2.  (a)  rectangular  cross  section,  (b)  to  ( j)  stress  strain  distributions  across  the  cross 
section  for  different  cases. 


Closed-form  Moment  Curvature  Response 


The  development  of  stress  strain  across  the  section  by  increasing  the  normalized 
compressive  strain  is  presented  in  Figure  3.  Stress  strain  develops  at  least  to  stage  4 
where  compressive  and  tensile  failure  is  possible  if  Xmax  =  jUuc  in  case  6,  or  Amax  =  F  in 
case  4.  Moving  through  different  stages  in  Figure  3  depends  on  the  controlling  value  for 
/.max-  Using  the  auxiliary  points  defined  in  Table  9,  the  transition  points  defined  as  tpn 
between  different  stages  in  Figure  3  could  be  presented  by  the  following  equations. 


Stage  4 

A.max=|rUc  >H|  Failure 


8C 


Failure 


Figure  3.  Strain  development  in  the  cross  section  at  the  different  stages  of  loading 


tpl2=Min(juc0,A ) 

tp2 3  ,  C)  orMin(/Jci ,  B ) 
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Where  indexes  i  and  j  refer  to  origin  and  destination  stages  respectively.  The  net  force  is 
obtained  as  the  difference  between  the  tension  and  compression  forces  in  Table  3, 
equated  to  zero  for  internal  equilibrium,  and  solved  for  the  neutral  axis  depth  ratio 
defined  as  k.  The  expressions  of  net  force  in  some  stages  result  in  more  than  one 
solutions  for  k.  Using  a  large  scale  of  numerical  tests  covering  possible  ranges  of 
material  parameters,  the  solution  of  k  which  yields  the  valid  value  0  <  k  <  1  was 
determined  and  presented  in  Table  10.  Moment  expressions  are  obtained  by  taking  the 
first  moment  of  the  compression  and  tension  forces  about  the  neutral  axis.  Curvature  is 
calculated  by  dividing  the  top  compressive  strain  by  the  depth  of  the  neutral  axis  tdi.  The 
closed  form  solutions  for  normalized  moment  M,  and  curvature  (pi  with  respect  to  the 
values  at  the  tensile  PEL  points  are  presented  in  Equations  5,  6  and  7  and  Table  10. 
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Where  Mpa  and  (Ppel  are  moment  and  curvature  (for  a  material  with  the  same  modulus  of 
elasticity  in  tension  and  compression)  at  the  tensile  PEL  and  are  defined  as 


M 


PEL 


bh  Es PEL 

6 


Ppel 


^ZpEL 

h 


(8) 


The  ultimate  moment  Mu  (flexural  strength)  is  computed  based  on  Equations  5,  8  and 
Table  10.  For  a  ductile  material  like  resin,  Mu  could  approach  Mx  at  very  large  X 
values.  The  normalized  moment  at  the  very  large  X  values  M  \  is  computed  by 
substituting  X  =  coin  the  expression  for  k  in  case  nine  of  Table  10  and  by  substitution  of 
X  =  oo  and  Kt  in  the  normalized  moment  expression.  Equations  9,  10,  and  11  present  the 

values  of  the  neutral  axis  depth,  normalized  moment,  and  curvature  for  very  large  X 
values. 

co 
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The  neutral  axis  depth  and  normalized  moment  are  a  function  of  characteristic  points  of 
tension  and  compression  stress  strain  model  (co,  y,p,pco,  jlici  )•  The  k9  and  M’g  expressions 
are  functions  of  all  characteristic  points  in  the  stress  strain  model  but  when  compressive 
strain  values  are  large,  they  are  independent  of  tensile  post  PEL  stiffness  (a)  and  strain  of 
ultimate  tensile  strength  (pti)-  For  an  elastic  perfectly  plastic  materials  with  equal  tensile 
and  compression  elastic  moduli  and  equal  yield  and  post  peak  flow  stress,  Equation  9  and 
10  yields  to  0.5  and  1.5  respectively,  validating  the  theoretical  value  that  the  plastic 
moment  capacity  of  a  rectangular  section  is  1.5  times  its  elastic  yield  strength  at  pco  =  Pci 
=  (Liti  =  1  presented  by  Salmon8  (1990).  For  a  set  of  parameters  y,p,pco,  Pci,  the  critical 
value  of  co  which  results  in  a  flexural  capacity  in  the  infinity  (before  failure)  greater  than 
the  flexural  capacity  at  the  tensile  PEL  point  can  be  found.  By  equating  the  normalized 
moment  for  large  top  compressive  strain  to  one  (M'x  =1 )  the  critical  value  of  post  peak 
tension  flow,  coCriticai,  is  expressed  as 

WcO  -  Pi1  c9  +  M  l 

3(///c0  -pfic0+pLicl)-\  (i2) 


Parametric  Study 

Parametric  studies  determine  the  sensitivity  of  the  different  parts  of  the  model  to  the 
relative  behavior  of  the  variables.  Although  polymer  materials  show  strain  softening 
behavior  with  a  percentage  of  the  ultimate  tensile  strength9,  a  complete  set  of  parametric 
studies  was  conducted  to  examine  the  strain  softening  materials.  Effect  of  the  tension  and 
compression  stress  strain  relation  on  the  response  of  polymer  materials  was  addressed  by 
studying  the  effect  of  variation  of  a  parameter  in  the  constitutive  law  while  all  other 
parameters  were  held  constant  to  the  typical  values.  Epon  E  862  epoxy  resin  was  chosen 
and  typical  mechanical  characteristic  of  E  862  were  extracted  from  Littell  et  al.6  (2008) 
and  are  shown  in  Table  11. 


Table  11:  Typical  characteristic  va 
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35 

The  flexural  strength  and  ductility  for  each  material  parameter  was  expressed  as  the 
normalized  moment  curvature  response,  independent  of  geometry  and  tensile  PEL 
strength.  Figure  4  shows  the  results  of  the  parametric  study  for  strain  softening  material. 
Figure  4  (al)  to  (a3)  presents  the  effect  of  constant  flow  tensile  strength  on  the  moment 
curvature  and  on  the  location  of  neutral  axis  depth.  co=  0.017  and  1.63  correspond  to 


constant  tensile  plastic  flow  equal  to  1%  and  100%  of  the  ultimate  tensile  strength 
respectively.  Figure  4  (a2)  depicts  that  moment  curvature  response  is  extremely  sensitive 
to  the  variations  in  constant  tensile  flow  as  the  location  of  maximum  flexure  and  the  post 
peak  regime  completely  changes  with®.  The  flexural  response  changes  from  brittle  to 
ductile  behavior  as  co  increases  from  0.017  to  1.63.  For  the  parameters  given,  Equation  12 
yields  to  oocriticai  =  0.39.  Figure  4  (a2)  also  shows  that  in  order  to  obtain  the  bending 
moment  at  large  top  compressive  strains  equal  to  or  greater  than  the  bending  capacity  at 
the  tensile  PEL,  the  required  tensile  plastic  flow  should  be  equal  or  greater  than  co=  0.39 
(25%  of  the  ultimate  tensile  strength).  When  co=1.42  the  response  characterizes  exactly 
the  material  behavior  of  Epon  E  862,  and  Equation  (10)  gives  M'x  =2.58  as  it  can  be 
seen  from  Figure  4(a2).  Figure  4  (a3)  shows  that  decreasing  the  level  of  tensile  flow 
decreases  the  neutral  axis  depth,  especially  for  co  values  less  than  0.5  (30%  of  UTS). 
Tensile  failure  was  the  governing  mechanism  in  all  the  material  responses.  Materials  with 
co  =  0.08  and  co  =  0.017  didn’t  experience  compression  yield  and  their  stress  strain 
relationship  always  were  in  the  elastic  compression  region.  This  is  the  reason  that  their 
neutral  axis  depth  and  moment  capacity  drops  sharply  by  applying  top  compressive 
strain,  and  they  show  a  brittle  material  behavior.  Table  9,  Equation  4  and  Figure  3  show 
that  the  stress  development  in  the  cross  section  is  totally  independent  of  tensile  plastic 
flow  in  stages  1,  2,  3.  Stage  4  would  also  be  independent  of  co  if  tp23  ^  C  in  stage  3.  By 
calculating  the  transition  points,  one  can  find  out  that  tpn  =  A  =  0.9167  in  case  2,  Ip 23  = 
juco  =  0.93  in  case  5,  tp34  =  E  =  3.358  in  case  7,  for  all  co  values.  For  co=  0.017  and  0.08, 
tp45  =  H  (//( co=0.01 7)  =  3.43,  and  H( co=0.08)  =  3.70)  and  tensile  failure  happens  in  stage 
5.  However,  for  co  =  0.393,  0.5,  0.83,  1.16,  1.42,  1.49,  and  1.62,  tp45  would  be  equal  to 
Pci  =  4.49  in  case  9  and  by  comparing  /(co)  with  pUc,  all  will  fail  in  tension  in  stage  6.  For 
co=  0.017  and  0.08  and  in  stage  5,  since  H  (co)  is  less  than  pcl,  tensile  failure  happens  and 
top  compressive  strain  is  less  than  pcl.  Figure  4  (bl)  to  (b3)  depicts  the  effect  of  UTS  on 
the  moment  curvature  and  neutral  axis  location  at  constant  tensile  post  PEL  slope  and 
constant  tensile  flow  stress.  Figure  4(b2)  reveals  that  an  increase  in  pti,  increases  flexural 
strength  and  this  is  clear  for  high  pn  values.  However,  the  amount  of  M'x  is  not  affected 

as  much  as  the  flexural  strength  since  for  p,j  >  6  the  moment  at  infinity  is  less  than  the 
flexural  strength.  Figure  4  (b3)  shows  that  by  increasing  the  ultimate  tensile  strength, 
neutral  axis  moves  downward  and  approaches  toward  k  =  0.5.  One  can  look  at  the 
variation  of  p,i  as  the  variation  of  CYS  to  UTS  ratio  which  is  defined  as 

YVco  +  Pit1  a — /fco)  ancj  ky  substituting  y=  1.19,  /?  =  0.3,  or  =0.24,  pc0  =  0.93,  pci  =  4.49, 
\+a\ptl  -1) 

it  is  clear  that  changes  in  ptj  from  2.75  to  8  will  change  the  CYS  to  UTS  ratio  from  1.53 
to  0.81.  Figure  4  (cl)  shows  the  compression  and  tension  model  with  the  post  tensile  PEL 
slope  and  strain  of  UTS  points  varied  from  0.15  to  0.6,  and  5.33  to  2.083  respectively,  at 
fixed  UTS  and  post  peak  tensile  flow.  Stress  strain  models  of  three  sets  of  a  and  pti  are 
shown  in  Figure  4  (cl).  Figure  4  (c2)  reveals  that  changes  in  parameters  a  and  p,i  slightly 
affect  the  moment  but  extremely  affect  position  of  the  neutral  axis  for  a  wide  range  of 
normalized  top  compressive  strains  between  1  and  4  and  that  will  change  the  stress 
distribution  across  the  section  between  elastic  and  post  peak  range.  In  order  to  study  the 
effect  of  post  compressive  PEL  stiffness,  the  range  of  parameters  /?  and  pcl  were  used  to 


represent  the  variation  in  (3  while  CYS  was  constant.  The  compressive  and  tensile  stress 
strain  models  are  shown  in  Figure  4  (dl).  Unlike  the  post  tensile  PEL  stiffness,  analyses 
of  Figure  4  (d2)  depicts  that  flexural  strength  is  quite  sensitive  to  the  variations  in 
parameter  [3  as  it  significantly  affects  the  flexural  strength.  Increasing  (3  and  decreasing 
jUd  slightly  affect  the  ductility  as  it  is  shown  in  Figure  4  (d2).  Figure  4  (d3)  shows  the 
profile  of  the  neutral  axis  position  versus  the  applied  top  compressive  strain.  Curves  of  (/ 3 
=  0.1,  jUd  =  11.61)  and  ( (3  =  0.2,  juci  =  6.27)  have  completely  different  shape  comparing 
to  others.  After  X  exceeds  the  compressive  PEL,  the  neutral  axis  depth  increases  sharply 
to  statically  equilibrate  the  axial  forces  in  the  cross  section  but  since  CYS  is  greater  than 
UTS  and  tensile  plastic  flow,  it  starts  to  decrease  in  post  peak  regions.  Material  with  (/3  = 
0.1,  jUci  =  11.61)  did  not  experience  any  yield  in  compression  and  all  are  failed  in  tension. 
Figure  4  (el)  shows  the  compressive  and  tensile  stress  strain  relationship  with  different  [! 
and  constant  pci.  CYS  is  from  62.04  MPa  for  p  =  0.1  to  137.54  MPa  for  P  =  0.6.  UTS  is 
70.2  MPa  and  constant  for  all  cases.  Increasing  the  compressive  yield  strength  100%,  the 
flexural  capacity  increases  around  40%  as  it  is  shown  by  the  normalized  moment 
curvature  plots  in  Figure  4  (e2).  Figure  4  (e3)  depicts  that  for  p  =  0.1  and  0.2  as  CYS  is  in 
the  range  of  UTS,  the  neutral  axis  is  almost  in  the  middle  of  the  section  along  all  loading 
stages.  Analyses  indicate  that  tension  is  the  governing  failure  mechanism  in  all  cases 
while  compression  strain  exceeds  the  yield  point.  In  order  to  study  the  effect  of  initial 
compressive  stiffness,  a  range  of  parameters,  y  and  pc0,  were  coupled  to  represent  the 
increase  in  relative  compressive  to  tensile  stiffness  from  0.8  to  1.4  at  a  fixed  compressive 
PEL  to  tensile  PEL  ratio  of  ypc()  =  1.107  as  shown  in  Figure  4  (fl).  Figure  4  (f2)  and  (f3) 
reveal  that  changes  in  the  relative  stiffness  slightly  affect  the  moment  curvature  and  the 
location  of  the  neutral  axis  especially  in  the  nonlinear  phase.  Analyses  show  that  all  the 
cases  experience  yielding  in  compression  and  failure  in  tension. 
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Figure  4:  Parametric  study  of  a  typical  strain  softening  material 


Load  Deflection  Relationship 

When  a  beam  is  loaded  beyond  the  modulus  of  rupture  (MOR)  in  a  material  with 
deflection  softening  behavior,  increase  the  deformation  decreases  load  in  the  distinct  zone 
(around  the  loading  nose  in  3PB)  in  the  cracking  region  while  the  rest  of  the  beam 
undergoes  unloading.  The  length  of  the  localized  zone  in  3PB  with  notch  or  groove  was 
assumed  to  be  1  mm  (around  2%  of  the  beam  length).  Moment  distribution  along  the 


Stress  (MPa) 


beam  is  obtained  through  the  static  equilibrium  of  the  beam,  and  the  curvature  of  each 
point  along  the  beam  is  obtained  though  the  moment  curvature  diagram.  For  epoxy 
resins,  if  compressive  modulus  of  elasticity  and  CYS  are  greater  than  the  tensile  modulus 
of  elasticity  and  UTS,  the  shape  of  the  moment  curvature  diagram  greatly  depends  on  the 
value  of  the  post  peak  tensile  stress.  In  order  to  obtain  the  load  deflection  response  for 
3PB  from  the  moment  curvature  diagram,  an  array  of  discrete  load  steps  is  defined  for  a 
given  moment  curvature  diagram  using  static  equilibrium.  The  specimen  is  loaded  from  0 
to  P,nax  in  the  ascending  portion  of  the  moment  curvature  diagram  from  0  to  Mmax.  The 
curvature  for  this  portion  is  determined  directly  from  the  moment  curvature  diagram. 
Once  the  moment  reaches  MOR,  the  curvature  distribution  along  the  beam  in  the 
softening  regime  depends  on  the  location  and  the  history  of  the  strain  at  that  point.  In  this 
study,  the  length  of  the  localized  zone  is  taken  as  1  mm,  thus  the  majority  of  the  sections 
along  the  beam  will  undergo  unloading.  For  the  sections  with  moment  less  than  the 
moment  at  limit  of  proportionality  (LOP),  the  curvature  unloads  elastically,  while  for  the 
sections  that  have  been  loaded  beyond  Mlop,  the  unloading  curvature  will  depend  on  the 
recovery  percentage  of  curvature,  which  can  be  taken  as  zero  in  the  monotonic 
displacement  control  3PB  test.  Knowing  the  curvature  values  along  the  beam  for  each 
load  step,  one  can  calculate  the  deflection  at  the  mid-span  for  each  load  step  using 
moment  area  method  numerically. 

Experimental  Results 

Tension,  compression,  and  3PB  bending  tests  were  conducted  on  two  hydraulically 
driven  materials  test  systems  (MTS  4411,  MTS  syntech  1/S),  and  an  electrical  desktop 
load  frame  (Dual  Column  Servo-all-Electric  Frame  Model  #:800LE)  as  shown  in  Figure  5 
at  room  temperature  and  at  low  speed.  An  interface  load  cell  (interface  model  SM-1000) 
was  used  to  measure  the  axial  load  in  all  the  tests.  Digital  image  correlation  technique 
(ARAMIS  4M)  was  used  to  study  the  strain  fields.  This  technique  recognizes  the  surface 
structure  of  the  object  to  be  measured  in  digital  camera  images  and  allocates  coordinates 
to  the  image  pixels.  Stress-strain  relationship  could  be  found  using  this  technique  as  it  is 
capable  of  capturing  loads  off  of  the  test  machine  and  showing  loads  in  the  output  file. 
Dog  bone  samples  with  a  gage  length  of  14  mm  and  a  rectangular  cross  section  of  3.18 
mm  x  3.43  mm  were  selected  to  conduct  the  monotonic  tensile  tests  at  the  rate  of  59 
pstr/sec  and  493  pstr/sec.  Small  cubic  samples  (4  mm  x  4  mm  x  4  mm)  were  tested  under 
monotonic  compression  at  the  rate  of  493  pstr/sec.  Small  beams  with  the  width  of  4  mm, 
thickness  of  10  mm,  and  length  of  60  mm  with  a  groove  in  the  middle  of  the  beam  were 
selected  to  conduct  3PB  tests.  Flexural  tests  were  done  with  different  speeds  of  0.217 
mm/min,  0.542  mm/min,  and  1.813  mm/min  approximately  corresponding  to  axial  strain 
rates  of  59  pstr/sec,  148  pstr/sec,  and  493  pstr/sec.  Epoxy  resin  Epon  E  863  with  a 
hardener  EPI-CURE  3290  using  a  100/27  weight  ratio  was  used.  Figure  6  illustrates  the 
experimental  true  stress  strain  curves  from  tension  and  compression  tests.  Through 
triggering  between  load  cell  and  optical  system  and  examining  the  strain  field  around  the 
loading  nose,  Moment  curvature  response  from  the  3PB  tests  were  obtained  and 
illustrated  in  Figure  7.  In  tension  and  compression  tests,  vertical  strain  along  the  samples 
was  calculated  by  taking  the  average  of  values  of  nine  stage  points  located  on  a  sign  “+” 


in  the  highest  strain  region.  Results  of  the  experiments  show  a  softening  behavior  beyond 
the  peak  point  followed  by  a  constant  plateau  before  failure. 

(a)  (b) 


Figure  5.  (a)  3PB  test  set-up;  and  (b)  compression  test 
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Figure  6.  (a)  True  stress  strain  tension  response;  and  (b)True  stress  strain  compression 
behavior 
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Figure  7.  (a)  Moment  Curvature  response  for  B_2m  samples  with  groove  from  3PB  test 
(b)  Load  deflection  response  of  B_2  and  B_2m  samples  with  notch  and  groove 


Flexural  Response  of  Epon  E  863 

A  simulation  was  used  to  study  the  moment  curvature  and  load  deflection  response  of 
Resin  Epon  863  and  to  evaluate  the  effects  of  out  of  plane  loading.  The  tension  and 
compression  stress  strain  response  and  the  simplified  model  for  strain  rate  of  493  pstr/sec 
are  shown  in  Figure  8  (a).  Figure  8  (a)  also  shows  the  modified  model  which  has  been 
obtained  from  back  calculation  of  the  experimental  moment  curvature  response.  The 
mechanical  properties  of  the  simplified  model  for  493  pstr/sec  are  presented  in  Table  12. 


Table  12:  simplified  stress  strain  mechanical  properties  of  Epon  E  863  at  493  pstr/sec 
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Figure  8:  (a)  Experiment,  model,  and  fitted  model  of  tension  and  compression  stress 
strain  curve  for  493  pstr/sec  (b)  simulation  of  moment  curvature  for  sample  B_2m-5  with 
groove  and  strain  rate  of  493  pstr/sec;  (c)  simulation  of  load  deflection  for  sample  B_2m- 
5  with  strain  rate  of  493  pstr/sec 
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Figure  9:  (a)  Experiment,  model,  and  fitted  model  of  tension  and  compression  stress 
strain  curve  for  493  pstr/sec  for  B_2-3;  and  (b)  simulation  of  load  deflection  for  sample 
B_2-3  with  strain  rate  of  493  pstr/sec 

Based  on  the  mechanical  properties  presented  in  Table  12,  the  nine  non-dimensional 
parameters  for  the  model  are  the  following:  pc0  =  1.148,  pci  =  3.53,  puc  =  15.70,  pti  = 
2.55,  put  =  21.54,  y  =  1.09,  a  =  0.395,  P  =  0.305,  and  co  =  1.376.  Figure  8  (b)  shows  the 
moment  curvature  curves  from  the  3PB  test  (sample  B_2m-5  with  groove)  compared  with 
the  simulation  results.  Figure  8  (b)  shows  that  the  tension  compression  model 
underpredicts  the  moment  curvature  response.  Two  main  reasons  for  the  under  prediction 
are  the  difference  between  stress  distribution  profiles  between  uniaxial  test  and  bending 
test,  and  the  effect  of  lateral  constraint  in  bending  tests.  In  tension  and  compression  tests, 
the  entire  volume  of  the  sample  is  subjected  to  the  same  load  and  has  the  same 
probability  of  failure.  However,  in  a  bending  test,  only  a  small  fraction  of  the  tension  and 
compression  regions  are  subjected  to  the  maximum  peak  stress.  Therefore,  the  probability 
of  crack  nucleation,  propagation,  and  failure  development  in  tension  and  compression 
samples  is  higher  than  the  bending  samples.  Table  13  compares  the  limit  of 
proportionality  (FOP),  modulus  of  rupture  (MOR)  of  samples  B-l,  B-2,  and  B-3  with 
PEF  stress  and  CYS  in  compression  of  sample  C-l  (a  representative  of  the  compression 
tests)  and  PEF  stress  and  UTS  in  tension  of  sample  T-l  (a  representative  of  the  tension 
tests). 


Table  10:  PEF,  CYS,  UTS,  FOP,  and  MOR  of  samp 
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Results  of  the  uniaxial  and  bending  tests  conducted  by  the  authors  showed  that  the  ratio 
of  FOP  to  MOR  is  up  to  14%  higher  than  the  ratios  of  PEFt  to  UTS  and  PEFC  to  CYS.  In 
bending  tests  without  axial  force,  while  a  part  of  the  cross  section  is  subjected  to 
compression,  the  other  part  is  subjected  to  tension.  This  stress  profile  results  in  an 
increase  in  tensile  and/or  compressive  modulus  of  elasticity.  To  quantify  these  effects,  the 
authors  propose  two  scaling  factors  to  modify  the  modulus  of  elasticity  and  strength  of 
the  material,  Emnd  =  C/  x  E,  and  £pEL,mod  =  £pel  x  C2/C1 ,  and  using  the  same  nine  non- 
dimensionalized  parameters.  Back  calculation  showed  that  Ci  for  Epon  E  863  is  around 
1.2  to  1.55  while  C2  is  between  1.05  and  1.15.  Various  amounts  of  imperfections  in 


material  directly  affect  C2  coefficient  but  maximum  of  ( 
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)  could  be  considered  as  the  upper  bound.  An  inverse  analysis  approach  of  flexural 


results  will  establish  a  statistical  relationship  between  the  compression,  tension  stress 


strain  and  the  flexural  response.  Inverse  analysis  showed  that  C/  and  C?  for  sample 
B_2m-5  are  1.2  and  1.12  respectively  as  shown  in  Figure  8  (b)  and  (c).  Figure  9  illustrate 
the  model  and  the  fitted  model  as  well  as  experimental  and  simulated  load  deflection 
response  for  sample  B_2-3.  The  rate  of  loading  in  the  3PB  test  for  these  sample  were 
1.813  mm/min  (for  B_2m-5)  and  1.76  mm/min  (for  B_2-3)  corresponding  to  493 
pstr/sec.  C/  =  1.55  and  C2  =  1.05  were  used  as  the  scaling  factors  for  sample  B_2-3  as  the 
result  of  inverse  analysis.  While  the  modified  models  in  samples  B_2-3  captured  the  load 
deflection  response  up  to  peak  load  well,  they  did  not  follow  the  deflection  softening 
behavior  up  to  failure.  Results  of  the  parametric  study  showed  that  simulation  of 
softening  can  be  improved  by  changing  the  post  peak  tensile  level  and  further 
adjustments  of  the  other  parameters.  Modification  of  tension  and  compression  models  in 
the  strain  softening  region  will  improve  the  simulation  of  the  model  in  the  post  peak 
response  for  bending  samples  with  notch. 

Concluding  Remarks 

Explicit  moment  curvature  equations  using  nonlinear  tension  and  compression  stress 
strain  relation  for  epoxy  resin  materials  have  been  developed.  A  multi-linear  stress  strain 
relation  for  epoxy  resin  materials,  consisting  of  a  constant  post  peak  response  in  tension 
and  constant  yield  stress  in  compression  has  been  used.  The  material  model  is  described 
by  two  intrinsic  material  parameters:  (a)  tensile  modulus  of  elasticity  and  (b)  tensile 
strain  at  the  PEL  point,  in  addition  to  five  non-dimensional  parameters  for  compression 
and  four  non-dimensional  parameters  for  tension.  A  parametric  study  showed  that  the 
normalized  moment-curvature  response  is  primarily  controlled  by  the  normalized  post 
peak  tensile  strength,  normalized  UTS,  normalized  post  compressive  PEL  stiffness,  and 
normalized  CYS.  Results  show  that  for  materials  with  small  values  of  post  peak  tensile 
strength,  reduction  of  moment  curvature  response  in  the  softening  regime  is  considerably 
fast  and  the  response  terminates  with  a  relatively  low  compressive  strain.  Materials  with 
higher  normalized  post  peak  tensile  strength  have  a  gradual  reduction  in  the  height  of  the 
compressive  zone,  therefore  larger  deformations  are  possible.  Results  show  that  while 
very  brittle  materials  have  a  moment  capacity  equal  to  or  less  than  the  moment  at  PEL 
point,  epoxy  resin  materials  with  a  considerable  amount  of  post  peak  tensile  strength  have 
a  moment  capacity  around  2.5  times  the  moment  at  the  PEL  point.  Results  showed  that 
for  resin  materials  with  post  compressive  PEL  stiffness  between  30%  and  40%  of  initial 
tensile  modulus  of  elasticity  with  a  fixed  CYS  strain  value,  the  moment  capacity  is 
around  2.5  times  the  moment  at  the  PEL  point.  Increasing  CYS  by  increasing  the  post 
compressive  PEL  stiffness  at  high  CYS  values  marginally  affects  the  moment  capacity  in 
polymer  materials.  Simulation  of  the  load  deflection  response  of  polymer  materials  in 
3PB  test  revealed  that  direct  use  of  tension  and  compression  data  under  predicts  the 
flexural  response.  Load  deflection  response  showed  that  the  nominal  flexural  stress  at  the 
peak  load  estimated  by  a  linear  elastic  mechanical  of  material  approach  is  higher  than  the 
uniaxial  ultimate  tensile  strength.  Curvature  distribution  along  the  beam  was  integrated 
up  to  the  mid-span  of  the  beam  to  obtain  load  deflection  response  in  3PB.  Simulations  of 
the  experimental  data  clearly  revealed  the  effect  of  stress  gradient  on  the  material 
behavior  as  the  uniaxial  tensile  and  compression  tests  yields  a  lower  tensile  and 
compression  strength  than  the  flexural  tests.  By  applying  two  scaling  factors  (Cj  and  C?) 


to  E  and  Spel,  tension  and  compression  stress  strain  model  presented  by  E,  £PEL  and  nine 
normalized  parameters  (juco,  juci,  juuc  ,  jUti  ,  juut  ,  7,  cc,  fi  and  co)  could  be  used  for 
simulations  under  flexural  loading. 
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Table  5:  Normalized  force  component  of  each  tension  sub-zone  for  each  case 


OJ  I  Is)  w|  to  w|  to 


Table  7.  Normalized  moment  arm  to  neutral  axis  of  each  tension  sub-zone  in  each  case 
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Table  8.  Normalized  moment  arm  to  neutral  axis  of  each  compression  sub-zone  in  each 
case 


Table  9:  Auxiliary  transition  points  in  the  stress  strain  development  diagram 
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Table  10:  Neutral  axis  depth  ratio,  and  moment  for  each  case 
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